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Abstract 

We present qualitative but precise description of the linear beam-size effect 
predicted for the processes in which unstable but long-living particles collide with 
each other. We derive physically pronounced equation for the events rate which 
proves that the linear beam-size effect corresponds to the scattering of one beam 
of particles on the decay products of the other. We compare this linear beam-size 
effect with the known logarithmic beam-size effect measured in the experiments on 
a single bremsstrahlung at VEPP-4 and HERA. 
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1 Introduction 



Let us consider a process with unstable but long-living particle(s) in the initial state. 
Calculating corresponding cross section, one can meet a singularity in the physical region 
related to the appearance of real (not virtual) particles in an intermediate state. The 
appearance of such singularities leads to a divergent cross section. To perform calculations 
in this case a careful analyses of the physical set-up is needed. 

Such situation takes place for a number of processes at muon colliders, for instance, 
— > eVeX, — > ev^X, ... which have been recently discussed in refs. [|I|, §. To 
be concrete, below we consider the process 

-> eu e W + . (1) 

One of the Feynman diagrams of this process which contains the abovementioned singu- 
larity (the so-called t-channel singularity in the physical region) is presented in the Fig.l. 
This diagram can be viewed as a sequence of two processes: 

H~ -> tVeV^ (2) 

and 

z/ At /i + - W + (3) 

both of which can occur for the real muonic neutrino. Therefore, there exists a region of the 
final phase space corresponding to the reaction (fj) where denominator of the propagator 
of the muonic neutrino q 2 can be equal to zero []. 




w 



Figure 1: The Feynman diagram for the reaction fi fi + — * ev e W + which has the t- 
channel singularity and gives the leading contribution in the region of small \q 2 \. 

In the paper [0] it has been shown that accounting for the finite sizes of the colliding 
beams gives a finite cross section for the processes with the t-channel singularity in the 
physical region. This result has been obtained on the base of the results of the paper || 
where a general technique for taking into account the finite sizes of the colliding beams 
has been worked out (for the review see H). The effective cross section obtained in 

dN I A\ 

da=— (4) 

1 For the appearance of such singularity it is necessary that both of the processes ([| ||) have a possibility 
to occur. 
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(where N is the number of events and L is the luminosity of the collisions) can be 

written in the form 

da = da st + da nst . (5) 

The cross section da st by definition corresponds to the wide region of the phase space 
defined by the inequality q 2 < — m 2 where m is the muon mass. This part of the cross 
section can be calculated by means of the standard rules of the relativistic scattering 
theory. For the discussed process the result of this calculation can be found in ||. The 
cross section da ns t corresponds to the narrow region of the final phase space where 

- m 2 < q 2 < t , t > (6) 

which includes the point q 2 = 0. It turns out || that the main contribution to this cross 
section is given by the pole of the propagator and that this contribution is proportional 
to the transverse size of the colliding beams a: 

da QSt oc a. (7) 

It is the linear beam-size effect. 

In the present paper we show that the non-standard piece of the cross section da QSt is 
defined by the sequence of processes (Q), ([$]) for real muonic neutrinos. Such interpretation 
of the results of ref. leads, in particular, to a more clear understanding of its region of 
applicability. 

In view of the realistic situation at future muon colliders, we assume that longitudinal 
I and transverse a sizes of the colliding beams satisfy the inequalities: 

E 

/< 7CT , 7 = — -, (8) 
mcr 

fl«cr (9) 

where r is the muon life time in its rest frame, cr = 660 m. Due to the condition (||]) 
the number of muons can be thought to be a constant during the time which is required 
for the colliding beams to cross each other. Note, that conditions (ffH© are perfectly 
fulfilled for the projects of muon colliders which are discussed in the literature 0. Using 
the beam parameters from ref. and the total energy 2E = 100 GeV one gets: 

a I 



10-*. (10) 

cr 7CT 

Subsequent part of the paper is organized as follows: in the next section we give a 
simple qualitative description of the linear beam-size effect; section 3 is devoted to the 
direct calculation of the number of events corresponding to the z/ M /i + collisions; in section 
4 we transform the result derived in ref. |2| for cr nst to the rest frame of the fi~ beam and 
compare it with the results of section 3; in section 5 we discuss a case when conditions 
(H)-(P) are not satisfied; our conclusions are presented in section 6. 



2 Qualitative description of the linear beam— size ef- 
fect 

We consider a collision of a single muon fi~ with the beam of muons of the positive charge 
H + . In this analysis it is convenient to use the rest frame of the /i - . Being unstable, 
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this muon is surrounded by a "cloud" of z/^'s which appear in the pT decay (certainly, 
there are also v e and e "clouds", but they are not interesting for us at the moment). The 
density of neutrinos in this cloud decreases as 1/r 2 with the growth of the distance r from 
the its time dependence corresponds to the exponential decay law and the angular 
distribution is isotropic. Therefore, one finds 

, , 6(t-r/c) ( t-r/c\ 
This density is normalized by the condition 

oo 

J cn„(r, t)Anr 2 dt = 1 
o 

which means that the total number of u^s which cross a sphere of radius r is equal to 
unity. 

We shall see below that the main contribution to the <r ns t comes from the distance 
r ~ a. The typical time of collision is therefore At ~ a/c. If the collision occurs at the 
time t which satisfies the inequality At <C t <C r, the neutrino density can be taken to be 
time independent in the collision region (these assumptions are justified by the obtained 
result): 

n„(r,t) =— !-^. (12) 

The distribution of the neutrinos over impact parameters g of the /i + pT collisions Q is 
given by 

d 2 g + f° dz 1 d 2 g 



-oo 



We note that the main contribution to this integral is given by the longitudinal distances 
of the order of 

\z\ ~ g. 

The number of muonic neutrinos which collide with the /i + beam of the radius a is 
therefore: 

f d 2 g tt a 
N v ~ — = --. (14) 

J ACTQ I CT 


The non-standard piece of the cross section of the process (|I|) is proportional to the 
number of neutrinos (|I4"D and can be estimated as: 

Onst ~ — (a u ^w) (15) 

CT 

where (<r Ufl ^w) corresponds to the cross section of the process (|3]) averaged over effective 
spectrum []. 

The estimate (|15D corresponds to the o" nst obtained in 0. The exact calculation is 
presented in the next section. 

2 The z-axis is antiparallel to the momentum p 2 of the fi + . 
3 The exact form of this spectrum can be found in ref. Qj. 
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3 Direct calculation of the number of events in z^/i + 
collisions 

Having performed the estimate, we proceed further with the exact calculation. We perform 
it in the rest frame of the \x~ beam in which neutrino density n v has the simple form ([12]). 
For simplicity, we neglect energy and angular spread of the particles in the muon beams. 
The number of events for the — > W + process can be written in the form 

dN v ^ w = dL vlx {u) a vll ^ w {s vli ). (16) 

Here dL Ufl (u) is the spectral luminosity of collisions. The cross section <j v ^^,^ 

depends on the quantity ^ 

Sup = (q + P2) 2 = 2uE 2 (1 + cos 6) + m 2 

where u> and 9, <p are the energy and the escape angles of z/ M , and E 2 is the energy of 

For unpolarized muons the energy distribution of the does not depend on the 
neutrino escape angles and has the well-known form B 



UJ 



2 ( 4uA du 



w{u)duj = l§[ — \ 3 — . (17) 

\mj \ m J m 

Let ni(ri) and U2(r2,t) be the particle densities of /i~ and /i + beams respectively. In 
accordance with eq. (|1^) the density of the neutrinos in the point r 2 is defined by the 
muons fi~ decaying at the point ri by the following equation: 

n 1 (r 1 )d 3 r 1 

dnJr 2 ) - 



47rcr(r 2 — ri) 2 

Therefore, the spectral luminosity of collisions is equal to (see, for example, |7j §12) 

dL v ^{uj) = w(uj)duj v rfny(r 2 )n 2 (r2, t) d 3 r 2 dt (19) 

where 

<ll>2 1 + cos^. (20) 



ujE 2 

Note that the vector r 2 — r\ is parallel to the neutrino momentum q: 

(r 2 -ri) || q. (21) 

In the next section we shall show that the result for dN QSt obtained in completely 
coincides with eqs. (|i6l) -([2~0"l). 



4 The number of events dN nst in the rest frame of the 
ji~ beam 

In this section we follow the line of reasoning of ref. [Q. The only difference is that we 
switch to the rest frame of the fi~ beam. It turns out that this reference frame is quite 
helpful for the physical interpretation of the results obtained in ref. |§. 

4 Below we use a system of units where h = 1 and c = 1. 
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We use the following notations (see Fig. 1): s = {p\ + p 2 ) 2 = 4E 2 is the square of the 
total energy in the center of mass frame, T = 1/r is the muon width, p\ = p\ = m 2 , P3+P4 
is the 4-momentum of the final e~v e system, y = (p 3 +p4) 2 /m 2 , q = pi —pz ~P4 = (u;, q) is 
the momentum transfer in the t-channel and x = qp 2 /p\p 2 . The square of the momentum 
transfer in the t-channel is equal to 

q 2 = -^—+t (22) 
1 — x 

where q_i_ is the component of the momentum q which is transverse to the momenta of 
the initial muons and to is the maximal value of q 2 

max{g 2 } = t = ^ - % - V) m 2 . (23) 

1 — x 

Note, that 

t > for y<l-x. (24) 

As is explained in the Introduction, usual approximations in the scattering theory 
(based on the collisions of the plane waves) result in the divergent cross section in this 
case. In contrast to this approximation, we use the formalism of the colliding wave packets 
(beams of particles) of macroscopic, but finite sizes which has been developed in ref. 0. 
In the overwhelming majority of experiments both descriptions provide the same result. 
However, in a number of cases the usual approximation is insufficient and consideration 
of the colliding wave packets is mandatory — see description of the experiments on a 
single bremsstrahlung at VEPP-4 || and at HERA || and the review H. The same can 
be stated about discussed processes with the singularities caused by real particles in the 
intermediate states (t-channel singularities in the physical region). The breakdown of the 
standard calculations clearly indicates inapplicability of the standard approximations in 
the scattering theory for these cases. 

Let us remind that in the standard approach the number of events N is the product 
of the cross section a and the luminosity L: 

dN = da L, da oc |M| 2 , L = v 12 J n x (r,t) n 2 (r,t)d 3 rdt (25) 

where V12 = |vi — v 2 | = 2 for the head-on collision of the ultra-relativistic beams. The 
quantities n,(r, t) are the particle densities of the beams. 

The transformation from the plane waves to the wave packets results in the following 
changes. The squared matrix element \M\ 2 with the initial state in the form of the plane 
waves with the momenta pi and P2 transforms to the product of the matrix elements Mfi 
and MJ V with different initial states: 

da oc \M\ 2 -> da(n) oc M fi M* fl , . (26) 

Here the initial state |z) is the direct product of the plane waves with the momenta 
ki = Pi + \k, and k 2 = p 2 — , while the initial state \i') is the direct product of 
the plane waves with the momenta k' x = pi — |k and k' 2 = p 2 + |k . Instead of the 
luminosity L the number of events starts to depend on the quantity 

L(r) = v 12 J m(r 2 - r, t) n 2 {r 2} t) d 3 r 2 dt (27) 
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through the following formula 



■^e»^)I(r). (28) 

Note that in this formula the densities n\ and n 2 are taken in different points - - it 
corresponds to the nonlocal interaction of colliding muons (contrary, in eq. (|25| ) these 
densities are taken in one and the same point). 

The characteristic values of k are of the order of the inverse beam sizes, i.e. 

1 

a 

Usually this quantity is much smaller than the typical scale for the variation of the matrix 
element with respect to the initial momenta. This is the case for dN st which is determined 
by the region of large values of —q 2 > m 2 . Here we can put k, = in da{n) which 
immediately results in the standard expression for the number of events (^5]). 

However, for the quantity dN QSt determined by the region (Q), the variable q 2 can go 
through zero. Therefore in this region q 2 changes considerably if initial momenta are 
varied on a quantity of the order oil/ a. Consequently, transformations (p6|)-(|2"8|) reduce 
to the following modification: 

1 i i 

(29) 



\q 2 + ie\ 2 t + ie t' — ie 
where 

q 2 = (pi-p 3 ) 2 , t = (h- P3 ) 2 , t' = (k[- P3 ) 2 . (30) 

Let us expand t and t' up to the terms linear in k. This gives 

t = q 2 -\(q 2 ), t' = q 2 + \{q 2 ) (31) 

where 

A(g 2 )=^Q, Q = -p 3 + ^Pi- (32) 
As the result, the divergent quantity 

to 



B 



dq 2 



\q 2 + ie\ 



transforms to 

r d 3 Kd 3 r ■ L(r) J dq 2 

J (2tt) 3 L J [g 2 - A(g 2 ) + ie] [q 2 + X(q 2 ) - ie] ' { ' 

—m? 

As has been shown in [0 the main contribution to the integral over q 2 in eq. ([33|) is given 
by the pole 

q 2 = -X(q 2 ) +iexi -A(g 2 = 0) + ie 

in the upper half plane. 
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In the rest frame of the /i , the vector Q is equal to 

Q = -p 3 = q = 

where u is the energy of the real (q 2 = 0) neutrino and n is the unit vector which defines 
direction of motion of the neutrino. 

Taking this pole and performing subsequent integrations over n and r, we obtain 

7r r L(rxi) , , , 

B = - -V^ dr. 34 
us J L 



o 



As a result 



dN nst = da vll ^w dT / L(m)dr . (35) 

uj J 
o 

For unpolarized muon beams 

r 

dY = -(1 -y)(l + 2y)dxdyd(p (36) 

7T 

where if is the azimuthal angle of the vector q. 
In the rest frame of the /i~ we have 

w =^(l-y), a: = -^- = i(l-y)(l + cosfl). (37) 

Here 7r — 9 is the angle between momentum q of the muonic neutrino and the momentum 
P2 of the positively charged muon. 

Using these equations we express differential width of the \i~ as 

dY = -!-(l -y) 2 (l + 2y)dydcos 
2nr 



Now adding all pieces together and using (see also eq. (|2lf) ) 

drd cos 9 dip = 
we finally get the following expression 

dN nst = g(y)dy J v 2^ n ^ v ^ f ) ^ r d 3 r 2 dt, 

g(y)= 2(1 -y) 2 (l + 2y). (38) 
It is easy to see that after redefinition 

2uj 

y = 1 , r 2 - r = ri 

m 

expression d38|) completely coincides with eqs. (p!6|)-(p0|). 
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5 Remark on the case of a large— size beam 



What happens if conditions (H)-(H) are not satisfied? Though this question is purely 
academic one for muon colliders, it can have a certain principal interest. There is no 
general answer in this situation and every such case should be considered separately. In 
this situation the method of preparation of colliding beams as well as decrease in the 
numbers of muons and neutrinos during the collision are of great importance. The true 
consideration in this case would require the modelling of how really unstable particle 
travels in the collider — it is clear that some characteristics of the collider are likely to 
enter the final result at this point. 

In this section we consider a particular case related to the paper [|lj . In this paper the 
divergent cross section of the process ([!]) has been regularized by introducing a complex 
mass for initial unstable particles: 

m — > m 1 . 

2 

Solving the energy-momentum conservation conditions it is possible to show then that 
the square of the momentum transfer in the t-channel gets small imaginary part: 

q 2 -> q 2 - l -m{l-y)T. 

Therefore, the singularity is regularized and the calculation of the cross section starts to 
be possible. Below we explain the physical meaning of the result obtained in this way in 
ref. 5 

Let us assume that \x beam has been prepared much earlier than the collision takes 
place. This means that almost all muons have managed to decay. In this case every \l~ 
produces a single with definite distribution over the energy fllTD and with the isotropic 
angular distribution in the rest frame of the \x~ . To obtain the results of ref. [I] we also 
have to assume that the fi + beam is so large that all neutrinos produced in the yT decay 
are still inside this beam having a possibility to collide with the muons of the positive 
charge. This situation is similar to the fast decay of the yT in a medium of the muons of 
the positive charge. In such situation the number of events is defined by the cross section 
<j m which corresponds to the cross section of the process @ averaged over the spectrum 
of neutrinos and their escape angles: 



smOdOdtp 

v w{uj)duj — (T V n-*w{Suit) (39) 



where v is defined in fl20|). Using eq. (]37D and integrating over y and ip we obtain 

i 

Ax(l — x){2 — x) a v ^w( xs ) dx. (40) 
This result is identical to the one which has been presented in ref. ffl. 



o 



3 See also ref. jnj, where a similar regularization has been used for the process nN* — ► ttN*. 
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6 Conclusions 



1. Comparing the results of sections 3 and 4, we conclude that the linear beam-size 
effect obtained in |2| corresponds to the scattering of the /x + beam on the "cloud" of real 
neutrinos -produced in the [i~ — > ev e v^ decay. Such interpretation assumes that the 
colliding beams have finite sizes and that decays of the negative muons on their way to 
the interaction point produce a cloud which fills the region occupied by fi + beam. 

Really amazing is the fact that within perturbative calculations these effects can be 
found only if one considers a collision of two wave packets which correspond to the beams 
of the colliding particles. Standard approximations in the scattering theory do not work 
in this case. 

It is interesting to estimate the "formation length" If of this cloud. In the center 
of mass frame the angular spread of the muonic neutrinos $ is of the order of $ ~ I/7. 
Hence one obtains 

a 

If ~ — ~ 7 a. (41) 
v 

For the discussed parameters of the beams f| and the total energy 100 GeV one finds 
If ~ 1 cm. 

From this we obtain the number of the neutrinos z/ M which participate in the collision. 
It can be estimated as 

N, -t- ~ N, — ~ 10 4 (42) 
C7 M r cr 

for Np ~ 10 12 . This number shows the level of statistical fluctuations which one expects 
for the collisions. 

2. It is interesting to compare the process (P which has been discussed in this paper 
and the processes with large impact parameters considered in the review || , for example, 
the process 

ep — > epj. (43) 
It is not difficult to show that the impact parameters g, which give essential contribution 



to the cross section of the process ([43]) in the standard calculations, are of the order of 



< 1 AE e E p (E e — E~ 



m e m e m p E^ 



Note that the quantity g m can be macroscopically large: for instance, for the collider 
HERA g m ~ 1 cm already for E 1 ~ 0.3 GeV. In this case there is a cloud of equivalent 
photons (EP) around every proton, these EP scatter on the electrons of the opposite beam. 
As a result, the interaction of electrons with protons becomes macroscopically nonlocal. 
In the discussed process ([!]) similar macroscopically nonlocal interaction between fi + and 
pT is due to the cloud of real u^s which appears as a result of the muon decay. Certainly, 
the nature of these clouds is quite different, however, one can describe both processes by 
means of the same technique. 

Distributions of EP and u^s in g are absolutely different for both of these processes. 
The cloud in the process (0) is distributed in the impact parameter space according to 



eq. (|T3|) , i.e. 

dN u oc d X 
Q 



10 



This distribution leads (after integration over g in the beam region, g ~ a) to the linear 
beam-size effect: 

da m oc a. 

The distribution of EP has another form 

dN EP oc ^ (44) 
g 2 

from which the known logarithmic beam-size effect immediately follows: 

da ep oc In a. (45) 

Let us emphasize that the logarithmic beam-size effect has been already established in 
experiments [^. 0. 

3. Let us finally compare the contribution of real and virtual u^s to the total cross 
section. First, we remind that because of the fixing of the region of the final phase space, 
the cross section cr st (see eq. (^) and the discussion after it) is completely determined by 
virtual z^'s with q 2 < —m 2 . Contrary, the contribution of real z/^'s in the intermediate 
state absolutely dominates in a nst - the relative contribution of virtual neutrinos to this 
piece of the cross section can be estimated as ~ 1/ (ay/to) ~ I /(am) ~ 10~ 10 . Let us 
stress, however, that the approach described in D allows to calculate the contribution of 
virtual neutrinos to cr nst as well as the contribution which appears due to the interference 
of real and virtual neutrinos. 
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